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PREFACE 


I 

This  Note  reports  on  a  portion  of  a  one-year  investigation  of 
strategic  communications  performed  for  the  Defense  Advanced  Research 
Projects  Agency  (DARPA1  under  RAND's  National  Defense  Research 
Institute,  a  Federally  Funded  Research  and  Development  Center  sponsored 
by  the  Office  of  the  Secretary  of  Defense.  It  presents  an  analysis  of 
the  distortion  experienced  by  a  hybrid  AM/PM  signal  that  has  been 
filtered  in  the  amplification  stages  of  a  typical  receiver  and  then 
demodulated.  The  findings  should  be  of  interest  to  designers  of 
communication  and  data  link  receivers  and  of  communication  satellite 
transponders.  The  research  was  performed  in  RAND's  Engineering  and 
Applied  Sciences  Department  as  an  element  of  the  Applied  Science  and 
Technology  Program. 
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I.  INTRODUCTION  AND  SUMMARY 
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The  AMBER  (AM  Broadcast  Emergency  Relay)  system  will  provide  voice 
and  data  communication  capability  throughout  the  continental  United 
States  for  use  during  crises  and  for  postattack  reconstitution. 

Existing  AM  broadcast  stations  will  be  pnase  modulated  to  form  a  highly 
connected  network  with  many  alternative  routes.  It  is  clear  that 
AMBER's  successful  implementation,  and  peacetime  or  crisis  operation, 
will  depend  critically  on  its  ability  to  operate  without  interfering 
with  conventional  broadcasts.  That  is,  except  in  certain  unusual 
circumstances,  listeners  using  conventional  AM  receivers  should  not  be 
aware  of  the  auxiliary  AMBER  phase  modulation. 

An  ideal  AM  receiver  will  pick  up  only  the  AM  component  of  a  hybrid 
AM/PM  signal  such  as  that  transmitted  by  an  AMBER  station. 

Unfortunately,  AM  receivers  typical  of  those  in  common  use,  even  if  they 
are  well  designed  and  properly  operated,  will  also  pick  up  the  PM 
component.  This  comes  about  because  the  selectivity  characteristics  of 
typical  AM  receivers  (as  determined  largely  by  their  intermediate 
frequency  amplifiers)  allow  phase  variations  of  the  input  signal  to  be 
transformed  into  amplitude  fluctuations  that  then  appear  in  the  receiver 
output.  This  is  true  even  if  the  IF  bandpass  characteristics  are 
symmetrical  and  the  receiver  is  properly  tuned  (though  these  conditions 
minimize  the  interference).  Only  in  the  case  of  an  ideal  receiver  with 
a  uniform  amplitude  response  and  a  linear  phase  characteristic  is  there 
no  interference. 

As  a  result,  some  interference  in  typical  AM  receivers  is 
inevitable  and  a  major  goal  for  AMBER  is  to  choose  a  data  rate  and  phase 
modulation  characteristics  such  that  the  resulting  interference  not  be 
perceptible,  even  to  a  critical  listener.  Experimental  and  theoretical 
work  to  date  suggests  that  signaling  rates  of  about  150  symbols  per 
second  can  be  accommodated  if  the  phase  modulation  waveforms  are 
properly  chosen. [1]  Advanced  modulation  schemes  capable  of  bandwidth 
compression,  i.e.,  sending  more  bits  per  symbol,  have  also  been 
investigated  [2]  to  improve  the  capabilities  of  AMBER.  The  general 
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approach  has  been  to  take  advantage  of  the  more  nearly  ideal  bandpass 
characteristics  of  typical  AM  receivers  near  the  center  of  their 
bandpasses.  Because  of  the  poor  audio  response  of  typical  AM  receivers 
(and  typical  listeners)  at  low  frequencies,  this  approach  allows  data 
rates  of  up  to  several  hundred  bits  per  second  without  interference. 

To  furnish  a  theoretical  basis  for  future  analyses  and 
calculations,  the  output  of  an  AM  receiver  to  which  the  input  is  a 
hybrid  AM/PM  signal  of  the  type  to  be  encountered  in  AMBER  is  determined 
and  presented  here.  The  system  to  be  analyzed  is  diagrammed  in  Fig.  ]. 
An  input  signal  is  assumed  to  be  simultaneously  amplitude  and  phase 
modulated.  The  spectral  response  of  the  receiver  is  represented  by  a 
bandpass  filter  that  may  be  unsymmetr ical  and  mistimed.  The  demodulator 
output  is  the  magnitude  of  the  filtered  input  for  an  AM  receiver,  or  its 
phase  for  a  PM  receiver.  (Although  the  principal  interest  for  AMBER  is 
in  AM  receivers,  the  results  for  the  output  of  a  PM  receiver  are  also 
presented  for  completeness.) 

The  results  are  in  the  form  of  an  expansion  in  which  the  leading 
term  is  the  linear  (i.e.,  undistorted)  reproduction  of  the  input  AM  and 
PM  and  the  successive  terms  represent  the  distortion.  In  both  cases, 
the  significant  effects  are  given  by  the  first  few  terms  of  the 
expansions  because  the  applications  of  interest  are  those  in  which  the 


Input 

signal 


Output 

signal 


Fig.  1 — Diagram  of  typical  AM  or  PM  receiver 


interference  is  small.1  Finally,  a  spectral  analysis  is  performed  on  the 
output  of  an  AM  receiver.  The  leading  DC,  linear,  and  nonlinear  terms 
of  the  output  power  spectrum  are  derived  by  assuming  that  the  amplitude 
and  phase  modulating  signals  are  uncorrelated  zero-mean  random  gaussian 
processes.  A  similar  analysis  for  a  phase  demodulator  remains  to  be 
done . 

This  work  is  an  extension  of  an  earlier  analysis  [3]  that 
considered  the  phase  demodulation  of  a  filtered  PM  signal.  The  results 
presented  here  are  generalized  to  the  use  of  a  hybrid  AM/PM  signal  and 
to  the  consideration  of  both  phase  and  amplitude  demodulation;  the 
current  approach  differs  in  that  it  uses  a  generalized  Taylor's  series 
rather  than  the  specialized  expansion  used  before.  Also,  as  an 
analytical  convenience,  it  uses  a  relatively  simple  univariate  expansion 
applied  to  an  exponential  form  of  AM.  The  results  are  then  reduced  to 
those  that  would  have  been  obtained  by  making  a  more  complicated 
bivariate  expansion  of  the  familiar  linear  form  of  AM.  Finally,  they 
are  specialized  to  the  case  when  the  receiver  is  properly  tuned  and  has 
a  symmetric  frequency  response. 

A  recent  extension  [4]  of  Ref.  3  is  related  to  the  work  presented 
here.  It  is  a  straightforward  extension  of  the  analytical  approach  used 
in  Ref.  3  to  the  case  of  a  hybrid  AM/PM  signal  input  to  an  off-tune 
receiver.  The  expansion  of  the  output  time  function  is  presented  for 
the  case  of  a  logarithmic  (rather  than  linear)  AM  envelope  demodulator 
and  a  linear  phase  detector.  No  spectral  results  are  presented. 

‘For  well  behaved  input  signals,  the  results  are  exact  if  the 
complete  expansions  are  used.  The  fact  that  each  expansion  is  a 
Taylor's  series  in  which  the  entire  linear  (or  undistorted)  component  of 
the  output  is  given  by  the  leading  term  suggests  that  the  distortion, 
which  is  given  by  the  remaining  terms,  is  well  approximated  by  the 
leading  nonlinear  term  when  the  distortion  is  small. 
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II.  THE  MODULATED  INPUT  SIGNAL 


When  written  as  a  real  function  of  time,  a  hybrid  AM/PM  signal  with 
a  linear  AM  component  becomes 


s(t)  -  [1  +  a(t)]cos[2*fut  +  *>(t)],  |  a(t  )  |  s  1 


(1) 


whe^.e  f^  is  the  carrier  frequency,  a(t)  is  the  amplitude  modulating 


function,  $(t)  s  the  phase  modulating  function,  and  the  restriction  on 
the  mdgnitur1  of  a(t)  prevents  overmodulation. 


If  a(t)  and  $(t)  are  slowly  varying  in  comparison  with  f  ,  the 


signal  is  said  to  be  narrowband  and  its  spectrum  S(f),  which  is  given  by 
the  Fourier  transform 


S(f)  -  f  s(t)e  ,2w"  dt 


(2) 


contains  most  of  its  energy  in  the  vicinity  of  +f  .  It  is  convenient 


for  analysis  to  consider  the  complex  signal 


ZUnV)  -  [1  +  a(t)]e‘*U)e‘2',J,  |a(t)|  s  1 


(3) 


which  is  approximately  analytic  and  has  negligible  spectral  content  at 
negative  f requencies . [ 5 ]  The  subscript  lin  denotes  that  the  signal 
corresponds  to  conventional  linear  AM. 


Generally  speaking,  if  the  spectrum  of  a(t)  is  band  limited  to  F  , 

3 


and  0(t)  is  constant,  the  spectrum  of  z, .  (f  )  is  limited  to  the  band  +F 

r  lin  -a 


about  f  .  However,  when  the  signal  is  phase  modulated,  its  spectrum 


extends  indefinitely  on  either  side  of  f  regardless  of  the  spectral 


j-  *_►  -■ 
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extent  of  $  because  of  the  transcendental  nonlinear  exponential 
operation.  Fortunate’y,  if  the  peak  phase  excursion  of  (»(t)  is  limited 
(to,  say,  +u/2  or  less)  and  the  spectrum  of  $(t)  is  band  limited  to  F  , 
the  spectrum  of  the  purely  phase  modulated  sigaal  (i.e.,  when  a(t)  is 
constant)  is  largely  confined  to  the  band  +F^  about  fQ.  If,  in 
addition,  the  waveform  of  $(t)  is  carefully  chosen,  the  spectrum  of  the 
modulated  signal  will  fall  off  rapidly  at  larger  separations  from  fQ. 
When  modulated  simultaneously  by  a(t)  and  $(t),  the  spectrum  is  extended 
somewhat  by  a  convolution  process  to  a  value  larger  than  that  produced 
by  either  phase  or  amplitude  modulation  alone.  Examples  of  this  are 
shown  in  Fig.  21  of  Ref.  1. 

A  hybrid  AM/PM  signal  using  exponential  AM  is  given  as  a  complex 
function  by 


_  e^e'2'Lt 


(4) 


where 


h(t)  =  a(t)  +  i<p(t) 


(5) 


and  where  the  subscript  exp  denotes  that  the  signal  corresponds  to 
exponential  AM.  Thus,  whereas  the  signal  given  by  (1)  is  a  function  of 
the  two  real  variables  a(t)  and  $(t),  the  signal  given  by  (4)  is  a 
function  of  the  single  complex  variable  h(t).  Inasmuch  as  the  expansion 
of  the  demodulator  output  is  obtained  using  a  Taylor's  series  (see  Sec. 
IV),  the  use  of  the  exponential  form  for  the  input  results  in  a 
univariate  expansion  in  terms  of  a  single  complex  variable  rather  than  a 
bivariate  expansion  in  terms  of  two  real  variables.  There  is  little 
difference  in  the  complexity  or  difficulty  of.  the  two  approaches  for  the 
first  term  or  two  in  the  expansion.  However,  higher-order  terms  in  the 
bivariate  expansion  grow  much  more  rapidly  in  complexity  than  those  from 
the  univariate  expansion.  In  fact,  general  recursive  relationships  can 


be  developed  in  the  univariate  case  to  permit  writing  terms  of  any 
desired  order  explicitly;  this  does  not  seem  to  be  possible  with  the 
bivariate  expansion. 

Of  course  the  analytical  simplicity  of  using  the  representation  for 
exponential  AM  is  of  little  value  if  the  results  so  obtained  cannot 
easily  be  related  to  those  for  linear  AM,  which  is  the  form  of 
modulation  of  practical  interest.  The  relationship  between  linear  and 
exponential  AM  can  be  seen  by  expanding  the  factor  exp  a(t)  in  (4)  in  a 
Taylor's  series  and  then  noting,  from  (3)  that 


Zexp(t)  =  [1  +  a(t )  +  0(a1)]ew~u)e‘~'1 1  =  z(m( f)  +  0(a2  )e,'eine 


,2 


2  2 

where  0(a  )  denotes  terms  of  order  a  and  higher.  Thus,  if  z  (t)  is 

exp 


used  as  the  input  signal  in  Fig.  1,  the  output  will  contain  the  desired 


result  represented  by  z^^Ct)  plus  extraneous  terms,  all  of  which  will 


involve  factors  0(a  ).  These  factors  are  readily  identified  and  deleted 
to  give  the  desired  result,  as  shown  in  Secs.  V,  VI,  and  VII. 
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III.  THE  FILTERED  AND  DEMODULATED  OUTPUT  SIGNAL 


Consider  an  arbitrary  linear  filter  having  an  impulse  response  g(t) 
and  a  steady-state  transfer  function  G(f).  These  are  then  a  Fourier 
pair 


gU)  -  f‘  G(f)eli'«  df 


(7) 


G(f)  =  J  g(t)e  ‘2'/l  dt 


(8) 


where  g(t)  is  real  and,  for  a  physically  realizable  filter,  vanishes 
when  t  <  0  (i.e.,  it  is  causal).  G(f)  is  complex,  in  general,  with  an 
even  real  part  and  an  odd  imaginary  part;  therefore 


G(f)  -  G  <-/) 


(9) 


where  the  overbar  denotes  the  complex  conjugate. 

The  filter  of  interest  here  is  a  narrow  bandpass  filter.  Such  a 
filter  is  called  bandpass  because  it  passes  a  band  of  frequencies 
preferentially;  i.e.,  the  magnitude  of  its  response  jG(f)|  is  large  in 
the  passband  and  small  elsewhere.  It  is  said  to  have  a  center  frequency 
f  ,  which  may  be  the  frequency  at  which  the  response  is  the  greatest,  or 
a  frequency  that  characterizes  the  center  of  the  passband  in  some  sense. 
The  choice  is  obvious  for  a  symmetric  filter  but  may  be  arbitrary  for 
others;  in  general,  it  suffices  that  a  center  frequency  be  specified, 
however  chosen.  If  the  width  of  the  passband  is  small  (by  any 


reasonable  measure)  in  comparison  with  the  center  frequency,  the  filter 
is  said  to  In'  narrowband. 


8 


Just  as  the  bandwidth  of  the  modulated  signal  is  related  to  the 
frequency  content  of  its  modulating  signal  relative  to  carrier 
frequency,  so  also  is  the  bandwidth  of  the  bandpass  filter  related  to 
the  apparent  frequency  content  of  its  impulse  response  relative  to  the 
center  frequency.  That  is,  the  impulse  response  can  be  written 


g(f)  =  r(t  )cos  [2xfct  +  0(0) 


(10 


where  r(t)  and  0(t)  are  real  functions  that  vary  slowly  in  comparison 
with  f  .  Furthermore,  just  as  the  spectrum  of  the  real  modulated 
narrowband  signal  contains  most  of  its  energy  in  the  vicinity  of  +f  ,  S' 
also  does  the  transfer  function  of  the  narrow  bandpass  filter  (having  a 
real  impulse  response)  contain  most  of  its  response  in  the  vicinity  of 
+fc-  Finally,  just  as  the  complex  narrowband  modulated  signal  z(t) 
given  by  (3)  is  approximately  analytic  and  has  negligible  spectral 
content  at  negative  frequencies,  so  also  would  the  complex  impulse 
response 


7 


( t )  =  r{t  )<’ 


tHU  I e‘-'l  I 


(11 


of  a  narrow  bandpass  filter  be  approximately  analytic  and  its  transfer 
function  T(f)  have  negligible  response  at  negative  frequencies.  (This 
suggests  the  possibility  of  an  analytic  filter,  which  is  described  in 
the  Appendix.)  It  should  be  noted  that  the  odd-conjugacy  relationship 
for  G(f)  given  by  (9)  is  not  valid  for  T(f).  Rather, 

\'(f  >  ~  0.  /  •  0  (12 

witli  the  bulk  of  the  response  of  Tlf  i  m  the  vicinity  of  f 


.  A.  -  *-  JW.,. 


When  the  complex  modulated  signal  given  by  (4)  for  a  hybrid  AM/PM 
signal  having  exponential  AM  is  applied  to  the  filter  having  the  complex 
impulse  response  y(t)  given  by  (11),  the  output  becomes 

x(t)  -  Jn'  ylrU(t  -  r)dr  =  £  r  (rh*1"1  '"ei2wLit  "dr 

=  e,2w,Jf'y  '(r)ehtt  "dr  (13 


where 


y’(t)  =  r(t)em,e,M/'  ,Jt 


and  (5)  has  been  used  to  introduce  h(t).  The  term  y’(t)  is  the  complex 
impulse  response  of  a  filter  that  is  the  lowpass  equivalent  of  the 
mistimed  narrow  bandpass  filter.  It  is  seen  from  (11)  and  (14)  that 


y  '(()  =  y(t)e 


from  which  it  follows  that  the  transfer  function  T'(f)  of  the  lowpass 
equivalent  to  the  mistuned  narrow  bandpass  filter  is  given  by 


n/)  -  X  y'(t)e  ,2nl,dt  =  Jo  y(t)e  ,MI,f  Hdt  =  T(f  +  /„) 


As  a  result. 


->  *> 


1.1 


7  V)  dt  =  r (/„) 


(17 


The  equivalent  lowpass  ri]ter  is  not  physically  realizable,  in  general, 
but  is  approximately  physically  realizable  for  a  nearly  symmetric1 
narrow  bandpass  filter  that  is  properly  tuned. 

When  the  filter  output  x(t)  given  by  (13)  is  amplitude  or  phase 
demodulated,  the  results  are,  respectively. 


f(t)  =  mag  x{t)  =  \/x (t )  x(t) 


(18 


or 


f(t)  =  arg  x(t )  =  Im  log  x(t) 


(19 


i2irf  t 

where  Im  denotes  the  imaginary  part  and  the  term  e  in  (13)  has 

been  suppressed.  For  convenience,  the  primes  in  (16)  and  (17)  will  be 
dropped  henceforth  and  the  quantities  Y(t)  and  T(f)  will  be  understood 
to  refer  to  the  impulse  response  and  transfer  function,  respectively,  of 
the  normalized  lowpass  equivalent  to  the  real  narrow  mistuned  bandpass 
filter. 


lTrue  symmetry  of  the  transfer  function  of  a  physically  realizable 
bandpass  filter  about  its  center  frequency  is  not  possible.  However, 
symmetry  can  be  well  approximated  if  the  bandwidth  is  sufficiently 
narrow. [ 6 ] 


.  r‘ -  r.  1'*  *  -  r. 


J  -r  s 


IV.  TAYLOR'S  SERIES  EXPANSION 


The  demodulator  outputs  given  by  (18)  and  (19)  are  vector  valued, 
i.e.,  each  input  modulation  function  h(t,u),  weft  defines  a  value  of  the 
output  that  is  a  point  or  vector  in  S2.  Such  vector  valued  functions  are 
treated  by  Hille  and  Phillips  [7]  who  show  that  under  certain  conditions 
the  Taylor's  series  expansion  of  a  vector  valued  function  f(y)  is  given 
by 


f{y  +  h)  =  Y-1  6nf{y;h) 

X,n[ 


(20 


where  6nf(y;h), 


the  nth  variation  of  f(y)  with  increment  h,  is  given  by 


6nf(y;h)  =  ^  f(y  +  ft) 


t-0 


(21 


It  is  also  shown  that  6nf(y;h)  is  a  homogeneous  polynomial  of  degree  n 
in  h  and  that,  in  particular,  6nf(y;h)  is  linear  in  h. 

Assuming  differentiability,  the  demodulator  outputs  given  by  (18) 
and  (19)  will  be  expanded  by  applying  the  above.  The  expansions  will  be 
about  the  point  y  =  0  and  h  will  be  taken  as  the  complex  modulation 
function  given  by  (5).  That  is,  f(y  +  h(t))  in  (20)  will  be  identified 
with  f(t)  in  (18)  and  (19)  by  setting  y  =  0  and  letting  h(t)  be  the 
variable  of  expansion  rather  than  t,  as  is  customarily  done  in  the 
Taylor's  series  expansion  of  a  time  function.  To  avoid  confusion  with 
exponents,  (20)  and  (21)  will  be  written  using  subscripts  rather  than 


V.  AMPLITUDE  DEMODULATION 


The  derivatives  of  the  amplitude  demodulator  output  (18)  are  best 
found  recursively  by  starting  with  the  square  of  (18).  Writing  £h  for 
t,  yields 


fHft)  =  x(ft)xUh) 


The  nth  derivative  of  the  left-hand  side  of  (24)  evaluated  at  4  =  0 
becomes 


dn 

dr 


f2m 


r-o 


rin 


fUh) 


£/«*> 

f-0 


r-o 


n 

2  nCr  A 

r-0 


where  (23)  was  used  to  introduce  6  and  C  is  the  binomial  coefficient 

n  n  r 


-tV  ‘ 


To  facilitate  differentiating  x(Ch),  let  (13)  be  written 


xtfh)  -  J  ye* 


where  (4)  has  been  used  for  z(t)  and  whore  a  convolution  integral  is 
understood.  The  nth  derivative  of  the  right-hand  side  of  (24)  then 
becomes 


dn 


x(th)x{th) 


dn 

r  yh  1  (h 

r  f  dn 

J  vn 

ye"  ye" 

>■ 
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where  the  subscripts  identify  the  variables  of  integration.  Equating 
(25)  and  (28)  then  yields 


n 

2  nCr  K  A 
r-0 


If  >  i  +  h2)n 


(29) 


which  is  the  basic  recursive  generating  formula  for  the  6  .  Evaluating 
(29)  for  successive  values  of  n  leads  to 

So  *  I  r()  | 

5, /So  ~  I\ 

S2/ 6o  =  1 2  -  (6]  /  5<,)2  =-  /,  -  /f 

S3/S0  -  h  -  3(5,  /«o)(«2  "So)  -  /:,  -  3/,  1 2  +  3/i' 

S4/S0  -  -  4(6,/  5()K63  /  6„)  3(6,  60)“  =  / 4  -  4/, I;,  +  1 8/ f / 2  —  3/ 2  —  15/f 

S5/S0  -  Is  -  5(5,  5o)(54  60)  -  10(5,  60)(S;,  5„) 

=  h  -  5/,/4  f  30 /,2/:,  -  150/ i'  /,  +  45/,/?  +  105/ ?  -  l0/2/;)  (30) 

and  so  on,  where  (17)  is  used  to  evaluate  6^  and  where 

In  !  ,  ('  f  m  -(h,  +  h2)n,  n  s  1  (31) 

2  I  I  0  I 

The  expansion  for  the  amplitude  demodulation  of  a  filtered  hybrid  AM/PM 

signal  having  exponential  AM  (4)  is  then  given  by  (22)  using  (30)  for 

the  6  and  (31)  for  the  /  . 
n  n 

REDUCTION  TO  LINEAR  AM 

As  mentioned  in  connection  with  (6),  the  results  obtained  using 
exponential  AM  can  be  reduced  to  those  for  linear  AM  by  deleting  terms 
0(a“")  in  the  expansion.  This  can  tie  done  by  noting  that 


15 


(h,  +  h2)n  -  Ha i  +  +  («2  -  «>2)1" 

-  ((a i  +  a2)  +  »(^>i  -  ¥>2)]"  -  l*(^i  -  *P‘i M n  +  n(°i  +  ®2)l*(*i  -*2)1" 

+  n{n  --  V(ai  +  a2)2[i(i^i  -  ^2)]n  1 

+  n-n - - —  (a i  +  a2)3[t(*»i  -  v>2)]n  3  +  '  •  ‘ 


(32) 


2 

It  is  apparent  that  terms  0(a  )  appear  in  all  but  the  first  two  terms. 
When  these  undesired  terms  are  deleted,  there  remain  only  the  two 
leading  terms  and  a  term  n(n  -  lla^a^  [ i ( 0 ^  *"  from  the  third 

term1;  the  fourth  and  succeeding  terms  vanish  in  their  entirety.  The 
result  is 


(h,  +  h2)n  -  [i(*i  -  *2)]n  +  ra(o,  +  a2)[i(*>i  +  "("  -  Daia2  [» (^>i  -  *>2)]"  2 

-  [n(ai  +  a2)  +  i(<pi  -  v,2)l l* («Pi  —  ¥’2)]"  1  +  n(n  _l)aia2[*(^i  _  ^2)]"  2  (33) 


Then,  (31)  becomes 


7172(^1  +  h2)  -  ^  ^  |  2  f  f  Trml(«i  +  a2)  +  ‘(vh  ~  ¥>2)]» 

7l72{[®(al  +  a2)  +  i(<fi  1  -  ¥’2)][,(¥’l  -  ¥>2)]n  1 

+  n(n  -  l)a,a2  [i(v>i  -  ¥>2)]"_2}. 


n  -  1 

(34) 
n  >  2 


lThe  term  containing  the  factor  a1  a2>  which  might  ordinarily  be 
2  2 

regarded  as  0(a  ),  is  not  0(a  )  in  this  context  because  it  results  from 
the  multiplication  of  terms  containing  the  factors  a and  a2>  which  are 

2 

to  be  retained.  The  terms  0(a  )  to  be  deleted  must  actually  contain 
2  3 

factors  a  ,  a  ,  etc.;  these  are  the  terms  deleted  in  (32). 


and  the  expansion  for  the  amplitude  demodulation  of  a  filtered  hybrid 
AM/PM  signal  having  linear  AM  (3)  is  given  by  (22)  using  (30)  for  the  6^ 
and  (34)  for  7  . 

SPECIALIZATION  OF  LINEAR  AM  TO  A  PROPERLY  TUNED,  NARROW, 
SYMMETRIC,  BANDPASS  FILTER 

In  this  case,  the  impulse  response  Y (t)  of  the  equivalent  lowpass 
filter  is  very  nearly  real  so  it  is  written  as  r(t)  to  accord  with  (11). 
Then,  noting  (17)  and  substituting  (34)  into  (30),  the  first  few  terms 
in  the  expansion  become 

=  J  T  =  To 

5,  =  /  ra 

h  =  (J  r<fi)  ~  f  r'f1 

53  =  3  [  —  /  ra*’2  +  f  r*  j ™ fra  ( f  r<>)  ]  (35) 

As  expected,  the  DC  term  consists  of  the  carrier  (transposed  to  zero 
frequency  by  (16))  passed  through  the  zero-frequency  response  of  the 
equivalent  lowpass  filter.  The  linear  term  consists  simply  of  the 
amplitude  modulation  filtered  by  the  equivalent  lowpass  filter.  The 
second-order  term  represents  distortion  produced  by  the  phase  modulation 
alone.  The  third-  (and  higher-)  order  terms  include  mixes  between  the. 
amplitude  and  phase  modulation.  If  there  were  no  phase  modulation 
(i.e.,  $(t)  =  0),  then  only  the  linear  term  would  remain. 


VI.  PHASE  DEMODULATION 


The  phase,  demodulator  output  (19)  is  found  most  conveniently  he 
expanding  log  x(t)  and  then  taking  the  imaginary  part.  Thus,  instead  of 
(22),  the  expansion  is  given  by 


f(h)  -  £  \  ^  *« 

„0  n 


(36) 


where  6^  is  again  given  by  (23). 


The  first  two  terms  of  (36)  can  be  obtained  directly.  Thus 


=  f  (fM 


log  x(fh) 


r-o 


f-o 


=  log  J  ye*  =  log  J  y  ^  log  T0  (37) 


r-o 


and 


a,  -  £  fm 


f-0 


d  f 


log  x(f h) 


|  y he* 


f=0  f  7P 


C  h 


f=0  Jt 


■  h  I* 


(38) 


where  (17)  has  been  used.  For  the  higher  derivatives,  it  is  better  to 
start  with  the  first  derivative 


^ ( Ch  1  lnM  X  ( (h  *  '  i-h 

a f  as  x ( s  h 


x({h ) 


\  .-.v, 


■V 


I 


and  write  it  as 
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Differentiating  n  -  1  times  more  then  yields 


n  - 1 


2 


(40) 


Now 


dn 

dr 


- 


dn 

dr 


i 


ye 


(h 


J  yh"  e** 


(41) 


so  using  (23)  to  introduce  6^  and  setting  £  =  0  permits  (40)  to  be 
written 


n  -1 

2  n-.Cj  yhrSnr  -  J  yh" 


(42) 


2 


„.,CrJ  yhr  6n_r-  f  yh" 


(43) 


where  the  first  term  has  been  extracted  from  the  sum.  Noting  (17)  and 
transposing  then  yields 


in  -  y;  [  f  yh"  ~  2  -1  Crin-r  /  yhr],  fl*  2 


(44) 


.s  . 


■  V  *.»  v  V  V  V  V  V  "J  '/  AV  v.v/r  W1T4  «r.  wv  rj  r. 


;■ 

E 


«  n  »t  *  ^  v  w  ^  ftvurva 


t  ■»  un  -  'amown  "■  »  ’  • 


“• 


S 

*- 


-  !3  - 


Letting 


«/n  *  ~  J  yh ",  »  >  1 


(45) 


the  leading  terms  in  the  expansion  become,  from  (37),  (38),  and  (44), 

So  -  log  r0 

Si-Ji 

5 2  “  e/  2  —  S\J  i  “  J 2  —  <7 '  \ 

63  —  t/3  —  2b?J  1  —  6je/ 2  “  e/3  —  3e/  ie/2  +  2J  j* 

64  •>  c/4  —  383c/ 1  —  362t/ 2  —  5je/ 3“  e/4  —  4e/  ie/3  +  12c/  jv/2  —  3e/|  —  6e/ ] 

85  “  e/5  —  \b^<J |  —  683c/ 2  —  452e/3  —  5je/ 4 

“e/5  —  5e/  je/  4  +  20e/  fe/  3  —  60e/  j*e/ 2  +  30e/(e/|  —  10e/-2e/3  +  24e/f  (46) 

and  so  on.  The  expansion  for  the  phase  demodulation  of  a  filtered 
hybrid  AM/PM  signal  having  exponential  AM  (4)  is  then  given  by  (36) 
using  (46)  for  the  6^  and  (45)  for  the  J^. 

REDUCTION  TO  LINEAR  AM 

As  was  done  with  amplitude  demodulation,  the  reduction  to  linear  AM 

o 

is  accomplished  by  deleting  terms  (Ha")  in  the  expansion.  Now 


hn  -  (a  +  i»"  -  («»"  +  na(i<fi)n  1  +  1  a2(iV)"~2  + 


(47) 


Deleting  terms  0(a  )  and  substituting  (47)  in  (45)  then  yields 


>.v.- 


v.  V.V 


/  ••  /  .  v  j- 


‘.O/?/  vV.v; 


V.'aV'. 


■LvLvNO 
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Jr,  -  Y~fy(na  +  *V)(*V)"  #il  (^8) 


Thus,  the  expansion  for  the  phase  demodulation  of  a  filtered  hybrid 
AM/PM  signal  having  linear  AM  (3)  is  given  by  (36)  using  (46)  for  the  6 
and  (48)  for  the  J^.  If  the  terms  in  a(t)  are  deleted  to  reduce  the 
results  to  the  case  of  pure  PM,  the  expansion  agrees  with  the  earlier 
one  given  in  Ref.  3.  The  first  few  terms  are 


Im  80  -  Im  log  T0 

Im  6:  -  Re  f 

l  o  J 

Im  &2-  ~  Im  1 

3 

Im  S3-  -  Re  [  jj-J 7v>3  -  A. f-w Jyp2  +  JL.  (f  w)  ]  (49) 

where  Re  denotes  the  real  part. 

SPECIALIZATION  OF  LINEAR  AM  TO  A  PROPERLY  TUNED,  NARROW, 
SYMMETRIC,  BANDPASS  FILTER 

In  this  case,  the  impulse  response  r(t)  of  the  equivalent  lowpass 
filter  is  very  nearly  real  so  it  is  written  as  r(t)  to  accord  with  (11). 
Then,  noting  (17)  and  substituting  (48)  into  (46),  the  first  few  terms 
in  the  expansion  become 


-  21  - 


Im  60  -  0 


Im  61  - 


Im  &2  - 


t ih 
tJ - 


Im  63 - J r<p3  -  [2 fra  fra*  -  fr^fr^2] 


To 


[3(jVa)  Jr^>  -  (Jrv>)  ] 


(50) 


As  expected,  the  DC  term  vanishes  because  the  zero-frequency  response 
of  the  equivalent  lowpass  filter  is  now  real.  The  linear  term  consists 
simply  of  the  phase  modulation  filtered  by  the  equivalent  lowpass 
filter.  The  distortion  terms  include  mixes  between  the  amplitude  and 
phase  modulation.  If  there  were  no  amplitude  modulation  (i.e.,  a(t)  = 
0),  the  second-order  term  Im  6^  (and,  probably,  all  higher-order  even 
terms)  would  disappear  but  the  third  order  term  Im  6^  (and,  probably, 
all  higher-order  odd  terms)  would  retain  intrinsic  PM  distortion.  Only 
in  the  case  of  a  white  filter  (infinitely  wide,  uniform-amplitude 
passband--in  effect,  an  all-pass  filter  but  with  a  possible  time  delay) 
does  the  distortion  disappear  completely.  This  can  be  seen  from  Im  6^ 
in  (50)  by  taking  r(t)  as  the  Dirac  delta  function,  which  is  the  impulse 
response  of  such  a  filter.  Then, 


Im  63  =  —  ^  +  ?>w~  —  2vw  “  0 


(51) 


as  indicated  above. 


/  /  A  /  .■  V  V  "  r  V  V  V  'jk  V  ’>  V  V  ">  .-  V%>. 
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VII.  SPECTRAL  ANALYSIS  OF  AMPLITUDE  DEMODULATOR  OUTPUT 


The  output  of  an  amplitude  demodulator  operating  on  a  filtered 
hybrid  AM/PM  signal  was  derived  in  Sec.  V.  The  leading  terms  in  the 
power  spectrum  of  that  output  are  developed  in  this  section  by  assuming 
that  the  amplitude  and  phase  modulation  components  (5)  of  the  basic 
hybrid  AM/PM  signal  are  uncorrelated,  zero-mean,  gaussian  random 
processes.  The  results  are  first  developed  for  the  case  of  the  complex 
exponential  AM  (4);  they  are  then  reduced  to  the  case  of  linear  AM  (3), 
and  finally  specialized  to  the  case  where  the  filter  is  symmetric  and 
properly  tuned.  The  expansion  for  the  amplitude  demodulator  output  is 
given  by  (22)  using  (30)  for  the  6^  and  (31)  for  the  I  ^ .  The  equivalent 
lowpass  filter  parameters  are  given  by  (15)  and  (16),  which  are  related 
according  to 


J  7 ,e*i2''‘dt  =  r„,  J  yte*iu*dt  -  Ttf 


(52) 


where,  as  before,  the  overbar  denotes  the  complex  conj’  •’ate  and 
subscripts  are  used  to  identify  the  variables. 

AUTOCORRELATION  FUNCTIONS  AND  POWER  SPECTRA 

The  power  spectrum  W(f)  of  a  function  f(t)  is  found  by  forming  the 
autocorrelation  function 


R  T  =  Eft  f (  +  r 


k53) 


where  E  denotes  the  expectation  operator,  and  then  taking  its  Fourier 
transforms  F 


W,  =  FR,  -  ['  R,e  ,2wl’  dr 


154 ) 


'  ' v‘>.\  •.  L-."  -.V.W  v' 


I'V.V.'J'  V.V.V.  /.  /.V  .•  V  -■  W  1-  .-f  -T  '  *  •  -  --  ' 
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where  R  and  W  are  a  Fourier  pair.  Let 


J,®6,  *  FE&M)6j(t  +  t) 


(55) 


where  the  subscripts  are  again  used  in  connection  with  the  6s  to 
identify  the  orders  of  the  terms.  Then,  from  (22), 


w/>-  2  Zw6,96 

i -0  ) -0  1 ' 


(56) 


which  is  the  desired  termwise  expansion  of  the  output  power  spectrum. 
It  should  be  noted  from  (55)  that  if 

5,®6,  "  Wijif ) 


(57) 


then 


ij9b,  =  FE6j(t) SM  +  r)  -  FEd,(t)6j(t  -  r)  -  FRi;(-t)  -  W^-f) 


(58) 


■w 

I 


In  forming  the  various  autocorrelation  functions,  it  will  be 
necessary  to  consider  the  expected  values  of  products  of  the  h  factors 
(5).  When  a  and  <t>  are  uncorrelated,  zero-mean,  gaussian  random 
processes,  it  follows  that 

Ea  =■  Eip  -  Ea<p  —  0 

Eh ,h2  -  E(ai  +  i*>,)  (a2  +  i* 2)  -  Ea,a2  -  E<p  m  “^21-  M-i  *  r2-j 

£/i,h2  -  £(a,  +  1*1)  (02  -  *>2)  -  ^ia2  +  £*W2  -  1  +  =  R 2+-i 

Ehxh2  -  £(a,  -  1V1)  (a2  +  ‘^2)  -  +  £*1*2  -  R°  1  +  =  R2+-i 

£h,7i2  -  £  (a ,  -  i<pi)  (a2  -  t>2)  -  ^aia2  -  £*>1*2  “  £2  -1  -  ££-1  -  £2-1  (59) 


% 

V 

v  v  v  v  V.'.'-V  v  ','.' 


where  (53)  has  been  used  to  introduce  R.  The  superscripts  in  (59) 
identify  the  functions  a  and  whose  autocorrelation  functions  are  being 
taken;  the  superscripts  +  and  -  refer  to  the  two  arithmetic  combinations 
of  these  autocorrelation  functions.  The  differences  expressed  as 
subscripts  refer  to  the  differences  in  the  arguments  of  the  variables  in 
question. 

To  consider  products  such  as  (h^  +  +  h^)n,  which  arise  when 

products  I  I  from  (31)  are  formed,  it  is  noted  that  for  an  even  number 
m  n 

of  random  variables 

Ef\f  i  ■  '  '  f  iN  ~  2  EfJjEfkf  i  '  ■  '  Efmfn  (60 

all 

pairs 


where  (i,j),  (k,l),  ....  (m,n)  are  N  pairs  of  integers  selected  from 
1,2,.,.,2N  and  the  summation  extends  over  all  the  1  35  (2N  -  1) 
possible  pairings.  When  the  number  of  random  variables  is  odd, 

Ef  if 2  '  '  '  f2N  i=0  (61 


In  applying  (60)  it  is  convenient  to  consider  +  h^  as  the  basic,  zero' 
mean,  gauss ian  random  process  and  to  note,  from  (59)  that 


E(h i  +  h2)(h,\  +  h 4)  *  R\  1  f  R\  \  +  R{  >  +  R 4  2 


(62 


For  convenience,  paralleling  the  notation  in  (59)  let 


Wf  =  W'(  •  Wf  FR, 


(63 


Additional  useful  formulations  obtained  from  (54)  are 


FRt.a.hRT.r.d  -  J  W^Wf  -*  +  «  -Mdp 


CROSS  TERMS 

When  the  cross  terms  6i®5j  are  formed  in  (56),  it  is  seen  from 
(30),  (31),  and  (55)  that  the  resulting  integrals  will  always  have 
products  of  hs  in  their  integrands  of  the  form 


(ft,  +  h2)a(h,  +  h4)b(hs  +  h6)c 


where 


a+6+c+  •••  =i4; 


Hence,  from  (61)  it  is  seen  that  all  cross  terms  for  which  i  +  j  is  odd 
will  vanish  when  the  expected  value  is  taken.  When  i  +  j  is  even,  the 
following  cross  terms  can  be  identified: 


i  +  j  0  :  6„®6„  FE  |  r„  | 2  =  |  r„  |  “6(0)  Leading  DC  term 


where  6(0)  is  t lie*  Dirac  delta  function 


i  +  j  -2  :  6„®6->  6»®60  Second-order  DC  terms 


Leading  linear  term 


(70 


5 1®6) 

t  +  j  -  4  :  60®54i  54®5u  Third-order  DC  terms  (71 ; 

5i®53,  63^6!  Second-order  linear  terms  (72 

52®52  Leading  nonlinear  term  (plus  others)  (73; 

The  DC  terms  (68),  (69),  (71),  and  portions  of  (73)  stem  from  pairings 
in  (66)  according  to  (60)  of  which  none  display  the  form  f  f  ^  The 

second-order  linear  terms  (72)  are  so  labeled  because  only  one  of  the 
pairings  is  of  the  form  f  f t_j_T *  thus,  they  represent  corrections  to  (70) 
which  is  the  leading  linear  term  consisting  of  a  single  f  f  ^  pairing. 
The  analytic  problem  thus  reduces  to  evaluating  6^05^  in  (70)  to  obtain 
the  leading  linear  term  and  the  appropriate  portions  of  6-0^2  in  Lo 

obtain  the  leading  nonlinear  terms. 

EXPONENTIAL  AM 
Leading  Linear  Term 

Using  (30)  for  6^  and  (31)  for  1  in  (70)  leads  to 
5,®«,  =  FE6$I  ,(t)1 2(t  +  r) 

-  FE - °  .  f  f  7u7.(h,  u  t  h,  ,) dude  |  )7u>,(h,  . ,  u  +  h, . ,  ,.)  du  dv  (74 

4 1  r„  1  J  J 

Expanding  and  applying  (62)  then  yields 
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5,8,6, - L  F  f  f  f  f  yuyvywyx(R ,  i+u  +  *+.,  +  R,-x*u)  du  dv  dw  dx  (75) 

4 1  r0 1  J  J 


The  Fourier  transform  of  the  autocorrelation  functions  is  found  by  using 
(64),  which  leads  to 


417^  M  y*™*y*W  ei2'f(-"'u)  +  Wf*ei2wfi~I+u) 


+  W/*ei2w,i  w*v)  +  Wfei2,/i-**u)]  du  dv  dw  dx  (76) 


The  transfer  function  of  the  equivalent  lowpass  filter  is  introduced  by 
applying  (52)  to  yield 


«1®«,  -  -~-j Tor, r0  +  w^r_,Tor0r^,  +  wf*  r0rfr,r0  +  w,  ur^or.,] 

4  I  Fo  I 


I  r0 1 
1 


[Wf+ 1  r0 1 2(  i  T/ 1 2  +  i  r_,  | 2)  +  21*7  Ker£rfr_fi 


4 1  r0 1 


(77) 


Finally,  the  result  is  expressed  in  terms  of  the  power  spectra  of  the 
amplitude  and  phase  modulating  signals  themselves  by  using  (63)  to 
obtain 


1_ 

4|f„| 


5i®6i  "  ; .  ^-7? 


Wf 

Wf 


r0r,  ±  r0r^  i  ■ 


(78) 


where  the  matrix-multiplication-like  notation 


.N  A  -  *  -  -  A  .  •  -  *  A  -  •  A  A,  *j.  *  . 


,  •*.  •*.  •  •' 


-V/. 

.•  >  *  •  .•  V  V  V  *.•  •  •  V  .* 
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Zl  I  r0r,  ±  r0r./ 1 2  =  wf  \  r0r,  +  Ter.,  1 2  +  wf  \  f0r/  -  r0r  ,  i 2  (79) 

wjr 


has  been  introduced  for  convenience  Equation  (78)  is  the  desired 
expression  for  the  leading  linear  term  in  the  output  power  spectrum 
(56). 

Leading  Nonlinear  Term 

The  derivation  of  the  leading  nonlinear  term  Parallels  that 

of  the  leading  linear  term  «£“l  above  except  that  many  tedious 
algebraic  manipulations  are  involved.  For  simplicity,  the  various  steps 
will  be  outlined  and  intermediate  results  indicated  without  detailing 
the  intervening  steps. 

Using  (30)  for  62  in  (73)  leads  to 


«2®«2  -  FE6$[I2(t)  -  /?U)][/2(t  +  r)  -  /?(*  +  t)] 

-  FEtflhU)  I2(t  +  r)  -  /2(t)/?(t  +  t)  -  /?(t)/2(f  +  r)  +  /?(<)/?(*  +  t)] 

-  T,  -  (T,  +  T3)  +  T< 


The  four  terms  in  (80)  have  been  identified,  respectively,  as  Tj,  T2> 

T3>  and  T^  for  convenience  and  will  be  evaluated  in  succession.  See 
(82),  (90),  and  (95)  for  the  resulting  formulations  for  T^  T2  +  T.^,  and 

V 

Nonlinear  Term  TJ .  Using  (31)  for  I2  in  (80)  leads  to 

Tj  -  FE — ^—rffyu^ht-u+fit-v^dudv  [fyuyu(kt  + ,  u  +  ht+r  „)2  du  dv  (81) 

4 1  r0  r 

The  expected  value  of  the  product  of  the  four  tern...  represented  by  the 
parentheses  in  the  integrands  leads  to  the  sum  of  three  products  of 
pairs  of  expected  values  according  to  (60).  It  is  seen  that  two  of 
these  are  identical  leading  to 


T,  -  — ^ — ?F  J'J'J'J'  Tu7i>7»7x[®(k«  u  +  K  v)2E(h,+r-u  +  /i|+,_i)2 


+  2E2(h, _u  +  /it  _„)(/»,  +T.u,  +  du  du  dw  dx 


(82) 


The  first  term  in  (82)  is  clearly  a  third-order  DC  term  because  neither 
expected  value  leads  to  a  pairing  of  the  form  f  f^.  Thus,  only  the 
second  term  need  be  considered.  Applying  (62)  to  (82)  then  yields 


7u7u7u,7xF(«r--u,+u  +  «,+-x+u  +  Kr+-„+u  +  R,--X+V)2du  dv  dw  dx  (83) 


Expanding  the  parenthesis  leads  to  sixteen  products  of  autocorrelation 
functions  whose  Fourier  transforms  can  be  found  by  using  (65).  After 
extensive  manipulation  and  the  use  of  (52)  to  introduce  the  transfer 
function  of  the  equivalent  lowpass  filter,  the  result  becomes 


tx  -  —i—y  f  { W7W/-_,2fie(r02r,r_,  +  r,r_piv_,r,_,) 

2  I  M)  |  J 

+  W'Wf.'AReYoW.fY'-fY,  +  r.^f,.^,) 

+  w;wf*.„[  i  r0|2(  |  T/ 12  +  |  r./ 1 2)  +  2 1  r.p | 2 1  T/_p | 2] }  dp  (84) 

Finally,  the  result  is  expressed  in  terms  of  the  power  spectra  of  the 
amplitude  and  phase  modulating  signals  themselves  by  using  (63)  to 
obtain,  after  further  manipulation, 


-\  /  j-s  -•  .•  vv  vv  vwv  w  v  /  vwv  .■ 


‘J*  'J- 


j  V  vVv 


I 


I 
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i  r 

w?  Wf^ 

T\  -  $  I 

2w;  Wf  _, 

2 1  r0  r  J 

• 

Wt  Wf^ 

Tor,  ;  r  ,r0  '  rf_,r,  :  r_,i>_,|2dp 


(85) 


i 

.* 


r 

r 


\ 

« 


\ 

■ 

} 


where  the  notation  is  an  extension  of  that  used  in  (79). 

Nonlinear  Terms  T ^  and  T y  It  can  be  seen  from  (57)  and  (58)  that 
and  in  (80)  are  a  mirror-frequency  pair.  Therefore,  will  be 
evaluated  first  and  found  by  writing  -f  for  f  in  T^.  Using  (31)  for 
7^  and  I  leads  to 


•  2 

T2  -  FE~~~-  g  J fyu7u(ht-u  +  ht-v)2du  du [jj-yu >„(/», +  T_U  +  ht+T_„)  du  dv]  (86) 


8  I  r0 1 


Using  (60)  in  the  same  fashion  as  with  (81)  leads  to 


I 


■ 

f 

r 


J 


T2 


g  |  p  '  '  *  J'7u‘Yu7u/7x7y7xl^'(^f-u  “h  ^ t-v )  E  (Rt  +  r -ui  +  +  x-x )(^f  +  r-y  +  ^h+r-z) 

+  2E(ht-u  +  +T _x )E (h(  _u  +  h( _„)(/lf  +  r~y  i  ^(+r-x)l  ^  ^  (87) 


As  in  (82),  only  the  second  term  need  be  considered  because  the  first 
term  is  a  third-order  DC  term.  Applying  (62)  to  (87)  then  yields 


4 1  r0 1 


T Ff  '  f 7u7v7i«7x7y7x(fl  T-W+U  +  ^T+-X  +  U  +  Rr-UI+V  +  ^r-x  +  «) 


(Rr-y  +  u  "f  R+-z+u  Ri *-y  +  u  +  Rr-z+v )  ^ 


(88) 


Expanding  the  product  of  the  parentheses  leads  to  sixteen  products  of 
autocorrelation  functions  whose  Fourier  transforms  can  be  found  by  using 
(65).  After  some  manipulation,  the  use  of  (52)  to  introduce  the 


>  *  •  *->*  ‘  *-  ’j*  .»  ‘ -T  *_•  .•  *.* 


/  .*  /  .* 
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transfer  function  of  the  equivalent  lowpass  filter,  and  invoking  the 
mirror- frequency  property  of  and  T^,  the  result  becomes 


T2  +  T3 


2 1  r. 


p/fe(r«?r„r/  pr  ,  +  2 1 r0 1 2r„r  pr,  prp_,  +  r03r_prp_,r,) 

I'o  l  J 


+  2 w„  Wf\  „Re(  1  r0 1 2r0rprp  ,r ^  +  r2 1  r,  j  2rpr  p  +  r2  1  r,_,  1 2rpr_p  +  1  r0 1 2r0r  pr,  „r/) 


+  w;  wf\ ,  1  r„  1 2[2 1  r„  1 2 1  r„  , 


I2  +  Re{ r„r  prp  ,r  ,  +  roT,r, _,!>)]}  dp 


(89) 


Finally,  the  result  is  expressed  in  terms  of  the  power  spectra  of  the 
amplitude  and  phase  modulating  signals  themselves  by  using  (63)  to 
obtain,  after  further  manipulation. 


i  r 

W°p  Wf., 

T2  +  T3  -  — ±-_ f  • 

2  W°p  Wf.p 

2  !  r0 1  J 

W*  Wf_p 

«e(r_,r0  :r0i7  \  r_pr,_,  :  rp_,rp)(ror_p  *  rpro)(r0r,-f  *  r,_pr0)dp  (90) 


where  the  notation  is  again  an  extension  of  that  used  in  (79). 
Nonlinear  Term  T Using  (31)  for  I ^  in  (80)  leads  to 


I 


2  2  -|2 

T4  -  FE - — t[//t«7.(*»  U  +ht  v)dudu]  [  ffyuyu(ht+T-u  +  ht+T.„)  du  dv]  (91) 

16  |  To  i 


Using  (60)  in  the  same  fashion  as  with  (81)  leads  to 
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T<  - 


is  i  r0 1 


fFj  '  •  '  J  7u7i 


[£(/it_u  +  /it -uM^t -10  +  /if  x)E(ht¥~-y 


+  ht  +  r-t  Hht  +r  -  a  +  ht  +T  .  (j) 


+  E  (h(  _  u  +  Aj -u)(Af  +  r-y  +  T  -  z)E  (H[  __  u  +  h(  g)(^t+r-a  ^  t  r  ■  fl)  (92) 

+  E(At_u  +  ht -u)(^i  +  r-a  +  At  +  T-g)E (A,  +  A<  j)(At(.t_y  +  A(  +  T_2)]  du  •  •  •  dff 

As  in  (82),  the  first  term  is  a  third-order  DC  term.  Also,  it  can  be 
seen  by  symmetry  of  the  variables  that  the  second  and  third  terms  are 
identical.  Applying  (62)  to  (92)  then  yields 


^4  *  Q  .  r  .  g  E  f  ‘  f  7u7vTrwTi7y7l7alfl^T  a  +  u  "F  Er-0+u  +  ^r-a  +  u  +  Ef-p  +  v) 

O  I  1  0  I 


(Rr-y  +  ui  +  Et-z  +  ul  +  Rr-y+j  +  fir_2+x  )  du  •  •  •  d/8 


(93) 


Expanding  the  product  of  the  parentheses  leads  to  sixteen  products  of 
auto-correlation  functions  whose  Fourier  transforms  can  be  found  by 
using  (65).  After  some  manipulation,  and  the  use  of  (52)  to  introduce 
the  transfer  function  of  the  equivalent  lowpass  filter,  the  result 
becomes 


T*  -  gj^|7 J l w wr-» 2Re rPr-P<r04r,  , rf  +  i  r0 1 4fp  f  f f  p) 

+  w  -  wf_p\Re  I  r0 1 2f02rpr  p(  i  r„  t  i 2  +  |  rf  p  i 2) 

+  w;  w,%  |  r„ 1 4(  i  rp  ,  \ 2  +  |  r,_p  1 2)(  |  rp  | 2  +  i  r  p  i 2)]  dP  (94) 


Finally,  the  result  is  expressed  in  terms  of  the  power  spectra  of  the 
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amplitude  and  phase  modulating  signals  themselves  by  using  (63)  to 
obtain,  after  further  manipulation. 


T< 


8 1  r0  I 


J 


wi 


Wf  , 
2W“  Wf_„ 

wr , 


w* 


I r„i\,  :  r0r, i 2 1 r0r^/  \  r0r, 2 d„ 


(95) 


where  the  notation  is  again  an  extension  of  that  used  in  (79).  The 
power  spectrum  of  the  leading  nonlinear  term  is  given,  as  indicated  in 
(80),  by  the  sum  of  (82),  (90),  and  (95). 

The  Complete  Result 

The  power  spectrum  of  the  output  of  an  AM  demodulator  when  the 
input  is  a  filtered  hybrid  AM/PM  signal  having  exponential  AM  is  given 
by  (56),  where  the  leading  DC  term  is  given  by  (68),  the  leading  linear 
term  by  (78),  and  the  leading  nonlinear  term  by  (80). 

REDUCTION  TO  LINEAR  AM 

The  reduction  to  linear  AM  for  the  basic  demodulator  output  was 

2 

achieved  in  Sec.  V  by  identifying  the  undesired  terms  0(a  )  in  the 

formula  for  I  in  (31)  and  deleting  them.  This  resulted  in  new 
n 

expressions  (34)  for  I  valid  for  linear  AM.  In  reducing  the  power 
spectral  results  to  linear  AM  it  is  more  convenient  to  identify  and 
delete  the  undesired  terms  in  the  spectra  than  to  rederive  the  results 
using  the  reduced  form  for  /  . 

Leading  DC  and  Linear  Terms 

Inasmuch  as  the  reduction  to  linear  AM  involves  the  suppression  of 
2 

terms  0(a  )  and  these  can  only  appear  for  7^,  n  >  2,  it  is  clear  that 
the  leading  DC  and  linear  terms  given  by  (68)  and  (70)  are  the  same  for 
both  exponential  and  linear  AM. 
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Leading  Nonlinear  Term 

Of  the  components  of  identified  in  (80),  it  is  seen  that  only 

T^,  ^2*  3n<*  ^3  ^nvo^ve  1 2 '  ^ °  reduce  these  terms  to  linear  AM,  the 

procedure  is  to  note  that 


(hi  +  h2)2  -  [(<*1  +  i<p  1)  +  (02  —  i<n)]2 

“  a?  +  2a, a 2  +  af  +  2i(ai  +  a2)(>pi  —  <p2)  ~  (^1  —  w)2 


(96) 


2  2  “2 
Thus,  the  terms  a^  and  a^  should  be  suppressed1  in  (h^  +  h^)  .  This  is 

done  by  writing 


h  -  +  h2)2  -  (a?  +  a22)] 


(97) 


in  the  spectral  analysis  rather  than  using  (31)  for  I 

Non] inear  Term  T^.  Using  (97)  for  / ^  in  (80)  leads  to 


Ti  -  F£-  r  4  //7u7.,[(ft«-u  +  ht-v)2  -  (a2-u  +  a2-J]  du  dv 


//  7u7u[(^  +  ,-u  +  /»t+r-u)2  -  (a2+  r-u  +  a2+r-i,)]  du  dv 
-  4|  r0|  2F/J/j7u7t'7u,7xF[(^-u  +  ht„)2(ht+ r-w  +  ht+r-x)2 

(^t-u  ^ t-U )  (@t  +  T  —  W  ^  ®f+T-x)  (@t-U  i"  -v)(ht +T~w  ^  h(  +  T-x) 


+  (a2-u  +  e2-u)(o2+r-u,  +  a2tr-i)]  du  dv  dw  dx 


(98) 


-  2  -  2 

The  leading  term  F.(h  +  h  )  (h  ,  +  h,  )  is  the  complete  term 

t-u  t-v  t+x-w  t+t-x  K 


:The  rationale  for  retaining  the  a^  a^  term  in  (96)  is  discussed  in 
connection  with  the  reduction  of  (32)  to  (33)  for  linear  AM. 


i 

/  -.W  V  CKV.*  -.v.v/'-o 


.lsv-v-  .'•^‘->7 


for  Tj  for  exponential  AM  in  (81).  Thus,  the  remaining  terms  in  (98) 
represent  the  correction  required  to  reduce  the  results  to  linear  AM. 

To  determine  the  form  of  these  correction  terms,  first  note  that 
applying  (60)  to  the  appropriate  terms  in  (98),  all  of  which  are 
squared,  consistently  results  in  the  production  of  pairs  of  terms  the 
first  of  which  always  represent  third-order  DC  contributions  (in  the 
same  fashion  as  in  (82)  and  (87)).  For  example,  the  expected  value  of 
the  first  correction  term  that  appears  in  (98)  becomes 

E(ht.u  +  ht-jW+r-w  =  E(ht-.u  +  h,-v)2Ea?+r  w  +  2E2(ht„u  +  ht.ll)at+T.w  (91 

Furthermore,  the  remaining  (i.e.,  the  second)  term  in  the  right-hand 
side  of  (99)  reduces  to 

E2(ht  _u  +  ht -v)cit  +r -w  =  E  [(fl< -t,  +  i<fit-u )  +  (flf -v  —  - v )]®t  +t -u 

-  E2(at  -u  +  at .  t.  )a,  4  r  _  w  (1< 


because,  from  (59),  Ea$  =  0.  Thus,  h  can  be  replaced  by  a  in  the 
correction  terms  in  (98).  Making  these  adjustments  and  applying  (62), 
but  with  the  superscript  a  rather  than  +  or  - ,  then  brings  the 
correction  term  in  (98)  to  the  form 
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ATl  =  “  271^  F  MM  yuyuyuyA2R°-w.uR°  «.+l  +  2 R^x+uRar-x+v  +  2R°_W+UR°„X,U 


(K 


+  2 Rar-W  +  VR°  x.v  +  (Rr-Wi u)2  +  (Rr  Itu)2  +  (Rx-w*v)2  +  (fi?-*tu)2]  du  dv  dw  dx 


It  is  apparent  that  when  (65)  is  used  to  take  the  Fourier 


transforms,  the  results  will  display  only  the  spectral  products 


Also,  inasmuch  as  the  autocorrelation  variables  represented  by  the 
subscripts  in  (101)  always  involve  one  duplication,  using  (52)  to 
introduce  the  transfer  function  of  the  equivalent  lowpass  filter  will 


Therefore , 


cause  every  correction  term  to  display  one  factor  or  T  . 
it  is  seen  that  the  reduction  of  (85)  to  linear  AM  can  be  achieved  by 
simply  deleting  all  terms  of  the  form  W ‘V*  T  F.T  Tt  that  involve  a 
factor  rQ  or  F^.  The  result  is 


Tx  “  iTrTP o  i  r"  /r"  +  r  »vi  -|2 


2W°  Wf 
Wf  Wf 


r0r/  •  r  ,r0  r,  ,\\  i'  r 


:V  <.  l 


dp  (102) 


Nonlinear  Terms  T ^  and  Tj.  Using  (31)  for  I ^  and  (97)  for  I  in 
(80)  leads  to 
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T2  ”  FEa\  r  |6  f  f7u7»[(^t  u  +  ht-S  -  (a?  u  +  a,2.*)]  du  du 

°  I  1  o  I 

[/J7u7w(/l(t,  -u  +  *fr  t,)dudv] 


8  I  r0  | 


P  .  4  J'TuTuTu/T*  Ty  T  z  ^[(^f  u  7  v)  w  +  t  r 


-  (O,2  u  +  Q(2  „)(/lMr  u-  +  *!.  r  iWlK  y  + 


l)(^I  H  >  ♦  T  Z  ) 

j]  du  dz  ( 103) 


The  leading  term  E(h  +  h  )2(h_,  +  h  )(h  ,  +  h  ,  )  is 

t-u  t-v  t+t -w  t+x-x  t+t-y  t+T-z 

the  complete  term  for  T ^  for  exponential  AM  in  (86).  Thus,  the 
remaining  terms  in  (103)  represent  the  correction  required  to  reduce  the 
results  to  linear  AM.  Using  ari  argument  similar  to  that  following  (98) 
and  recalling  that  T^  and  T^  had  mirror  symmetry  in  frequency,  it 
follows  that  the  correction  term  can  be  brought  into  the  form 


A(T,  +  T;i)  = 
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As  was  the  case  with  AT^  in  (101),  applying  (65)  to  take  the 

Fourier  transform  will  result  in  the  production  of  only  the  spectral 

products  W3W''  Also,  inasmuch  as  (104)  is  a  sixfold  integral  and  the 
P  1  ”  P 

autocorrelation  variables  represented  by  the  subscripts  in  (104)  always 

involve  only  three  variables,  it  follows  that  when  (52)  is  used  to 

introduce  the  transfer  function  of  the  equivalent  lowpass  filter,  every 

correction  term  will  display  three  factors  or  Therefore,  it  is 

seen  that  the  reduction  of  (90)  to  linear  AM  can  be  achieved  by  simply 

deleting  all  terms  of  the  form  W  av/l  P  I"  T  P  c  T  P  that  involve  exactly 

p  f-p  a  B  y  6  e  4  } 

three  factors  P  or  P^.  The  result  is 
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The  power  spectrum  of  the  leading  nonlinear  term  is  given  by  the  linear- 
AM  version  of  (80);  that  is 


S2®&2  -  T,  (T,  .  7’;,)  4  T4  CIO 


See  C102),  f 105)  ,  and  (95)  for  the  appropriate  formulations  for  T^ ,  T^  + 
T  ,  and  T. . 


The  Complete  Result 

The  power  spectrum  W(f)  of  the  output  of  an  AM  demodulator  when  the 
input  is  a  filtered  hybrid  AM/PM  signal  having  linear  AM  is  given  by 
(56),  where  the  leading  DC  term  is  given  by  (68),  the  leading  linear 
term  by  (78),  and  the  leading  nonlinear  term  by  (106).  These  are 
repeated  here  for  convenience 


"W-  2  ZjtttW 
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(56) 


where  the  significant  terms  are  given  by 
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Leading  linear  term 


(78) 


<5'2®<52  =  T\  -  {T2  +  T.\)  +  T4  Leading  nonlinear^  term  (106) 
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6(0)  in  (68)  is  the  Dirac  delta  function.  The  subscripts  are  the 
arguments  of  the  functions  to  which  they  are  appended.  W  and  W  are 
the  power  spectra  of  the  amplitude  and  phase  modulation  functions, 
respectively.  The  mat r ix-mu 1 1 ipl icat ion- 1  ike  notation  in  (105)  and  (95) 
denotes,  from  (79), 


jnl  (a  >  b)  A  (a  t  b)  t B(a  -  b) 
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A 
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(a  '  6)(c  d)  -  A(a  +  6)(c  +  d)  +  B(a  +  b)(c  -  d)  +  C’(a  -  6)(c  -  d) 


SPECIALIZATION  OF  LINEAR  AM  TO  A  PROPERLY  TUNED,  NARROW, 
SYMMETRIC,  BANDPASS  FILTER 

The  spectral  results  for  linear  AM  will  be  specialized  to  the  case 
when  the  narrowband  filter  is  nearly  symmetric  and  properly  tuned.  This 
case  is  of  interest  because  it  corresponds  to  the  situation  where  an 
ideal  AM  receiver  using  an  envelope  detector  is  properly  tuned  to  a 
hybrid  AM/PM  signal.  Also  it  is  a  case  in  which  nonlinear  terms 
simplify  enough  to  permit  some  interpretation.  The  specialization  is 
done  by  noting  that  the  impulse  response  of  the  equivalent  lowpass 
filter  is  very  nearly  real  for  this  case;  that  is  0(t)  ~  0  in  (11). 

Thus,  its  transfer  function  displays  the  conjugate  symmetry  property 

(9)  of  the  filter  itself  so  =  T  ^  and  is  real. 

The  leading  DC  term  (68)  is,  of  course,  not  affected.  The  leading 
linear  term  (78)  demonstrates  an  interesting  symmetry  even  before 
specialization  to  a  properly  tuned,  symmetric  filter.  It  is  seen  from 
(79)  that  ,  the  spectrum  of  the  amplitude  modulating  signal,  is 
multiplied  by  the  square  of  the  magnitude  of  the  sum  of  and  F 

A 

whereas  ,  the  spectrum  of  the  phase  modulating  signal,  is  multiplied 
by  the  square  of  the  magnitude  of  their  difference.  Thus,  the 
appearance  of  the  phase  modulated  signal  in  the  output  of  an  AM 
demodulation  is  related  to  the  lack  of  symmetry  in  the  filter.  In  the 
limit,  as  ,  the  phase  modulation  disappears  and  (78)  becomes 


5,®5,  -  Wf  |  r,  I  2 


(107) 


which  agrees  with  the  expansion  given  by  (35)  where  it  was  noted  that 
the  linear  term  6j  consists  of  the  amplitude  modulation  filtered  by  the 
equivalent  lowpass  filter.  The  result  (107)  is  valid  for  both 
exponential  and  linear  AM. 
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For  the  case  of  exponential  AM,  the  leading  nonlinear  term  (80) 
simplifies  somewhat  but  the  results  are  of  little  practical  interest  and 
are  not  shown  here.  However,  for  the  case  of  linear  AM  the 
simplification  of  the  leading  nonlinear  term  (106)  is  dramatic.  Letting 
=  f  £  and  Tq  be  real  in  (102),  (105),  and  (95)  leads  to 


I 

■  which  is  in  agreement  with  the  expansion  given  by  (35)  where  it  was 

noted  that  the  leading  nonlinear  term  6^  represents  distortion  produced 
by  the  phase  modulation  alone.  A  simple  interpretation  of  '108)  does 
^  not  appear  possible  but  the  form  of  a  difference  between  and  l*  F  j 

suggests  that  the  integration  on  p  is  a  conditioning  of  the  phase 
modulation  spectrum  over  the  filter  dissymetry  as  a  function  of  all 
k  possible  center  frequencies.  That  is,  the  wider  the  filter  is  with 

|  respect  to  the  bandwidth  of  the  phase  modulation,  and  the  flatter  its 

response,  the  less  the  value  of  (108)  will  be. 
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VIII.  CONCLUSIONS 


The  analysis  presented  in  this  Note  is  purely  theoretical.  The 
intention  was  to  develop  the  formulations  that  can  be  used  in  subsequent 
numerical  calculations.  Principal  interest  was  focused  on  the  output  of 
an  AM  receiver,  which  may  be  mistuned  and  whose  spectral  response 
characteristics  may  not  be  symmetrical,  when  the  input  is  a  conventional 
AM  signal  that  is  simultaneously  phase  modulated.  An  expansion  of  the 
output  time  function  is  given  in  Sec.  V  and  the  output  power  spectrum  in 
Sec.  VII.  An  expansion  of  the  output  time  function  from  a  PM  receiver 
is  given  in  Sec.  VI  for  completeness  but  the  corresponding  output  power 
spectrum  remains  to  be  derived. 

The  situation  of  practical  interest  is  the  one  in  which  the 
undistorted  output  (i.e.,  the  linear  component),  which  is  given  by  the 
leading  term,  is  large  in  comparison  with  the  distortion,  which  is 
conveyed  by  the  remaining  terms  in  the  output  time- function  and  power- 
spectrum  expansions.  Thus,  the  principal  contribution  to  the  distortion 
is  contained  in  the  leading  distortion  term  if  the  input  modulating 
signals  are  well  behaved  and  the  expansions  converge  rapidly.  The 
output  time-function  expansion  probably  converges  for  modulating  signals 
of  practical  interest.  The  situation  is  less  clear  with  respect  to 
output  power-spectrum  expansion  because  it  was  necessary  to  assume  that 
the  modulating  signals  were  independent  random  gauss ian  processes  to 
make  the  analysis  tractable.  Some  questions  concerning  convergence  are 
discussed  in  Ref.  3  but  the  matter  is  by  no  means  resolved.  Thus,  the 
expansion  must  be  used  with  care;  it  is  undoubtedly  valid  only  for  the 
case  of  small  distortion  and  should  be  verified  experimentally,  if 
possible. 
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Appendix 


THE  ANALYTIC  FILTER 


Consider  a  real  physically  realizable  filter  having  an  impulse 
response  g(t)  and  transfer  function  G ( f )  defined  by  (7)  and  (8), 
respectively.  An  equivalent  analytic  filter  can  be  defined  having  an 
impulse  response 


u(t)  =  l-  [£(£)  +  ig(t)] 


(A .  1 ) 


where  g(t)  is  the  Hilbert  transform  of  g(t).[5]  If  the  Fourier 
transform  F  of  g(t)  is  given 


Fg(t)-G(f)  (A. 2) 


then 


Fg(t)  =  -i  sgn  f  G(f)  (A. 3) 


whe-  ; 

f  /.  /  >  0 

sgn  f  =  I  0,  {  =  0 

1-1,  f  <  0  (A. 4) 

Thus,  the  analytic  impulse  response  u(t)  is  that  of  a  filter  having  a 
transfer  function 

(Gif  ),  /  >  0 

U{f)  =  (o,  /-()  (A. 5) 
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Inasmuch  as  the  analytic  filter  has  the  same  transfer  function  as  its 
equivalent  real  filter  for  positive  frequencies,  its  response  to 
analytic  signals,  which  contain  only  positive  frequencies,  is  identical 
to  that  of  the  real  filter.  That  is,  if  s(t)  is  the  input  signal,  then 

J0  £(r)[s(£  -  r)  +  is(t  +  r)]dr  =  ^  J  [g(r)  +  i^(  r))  [.s' (t  -  t)  +  isit  r  )]dr  (A.  6) 

which  can  be  verified  by  taking  the  Fourier  transform  of  both  sides. 

The  analytic  filter  is  not  the  same  as  the  filter  having  the 
complex  impulse  response  y(t)  given  by  (11)  because  the  real  and 
imaginary  parts  of  Y(t)  are  not  a  Hilbert  pair. [8]  That  is, 

r(t)sin[2irfct  +#(£)]  /  Hr  (t  )cos  [2irfct  +(/(£)]  (A.  7) 

where  H  denotes  the  Hilbert  transform.  However,  the  relationship  is 
approximately  true  for  a  narrow  bandpass  filter,  which  is  the  reason  why 
the  response  T(f)  of  the  complex  narrow  bandpass  filter  can  be  assumed 
to  be  negligible  for  negative  frequencies . [9] 

It  should  be  noted  that  the  analytic  filter  is  realizable  in  the 
sense  that  the  real  part  of  its  impulse  response  is  the  causal  impulse 
response  of  a  physically  realizable  filter.  The  imaginary  part  of  its 
impulse  response,  which  is  the  Hilbert  transform  of  the  real  causal 
impulse  response,  is  not  causal  in  general;  hence,  the  extended  lower 
limit  ori  the  right-hand  integral  in  (A. 6).  Using  the  analytic  filter 
may  be  a  nuisance  for  time-domain  analysis  because  of  the  need  to  find 
the  Hilbert  transform.  Its  use  in  the  frequency  domain,  on  the  other 
hand,  is  simple  because  it  merely  entails  deleting  the  transfer  function 
of  the  real  filter  for  negative  frequencies. 
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